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Abstract 

We study the cosmology of quintessence models in an extended theory of gravity in Lyra's ge- 
ometry. By analyzing the possible interactions between the quintessence scalar and the intrinsic 
displacement field in Lyra's geometry, we obtain the closed form solutions of the modified Fried- 
mann equations for four classes of quintessence models. Though the presence of the geometrical 
displacement field promises the possibility for the effective equation of state uj of the quintessence- 
displacement mixture crossing the cosmological constant boundary, the reliable quintessence sce- 
narios in Lyra's geometry with stable perturbation modes are still those in which —1 < uj < 1. 
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INTRODUCTION 



Recent cosmic observations l|-|l4| have indicated that our Universe is undergoing an 
accelerated expansion at the present epoch. The cause for such a cosmic acceleration is 
attributed to an unknown dominant energy component, dubbed dark energy, with negative 
pressure generating thus repulsive gravitational forces that counteract the attractive forces 
produced by radiation, baryons and the cold dark matter. However, the exact nature of 
dark energy is currently a significant part of the realm of speculations. Some believe that 
dark energy is the energy of the quantum vacuum, modelled by the cosmological constant 
A of general relativity. Interpreting dark energy as a cosmological constant means that the 
density of dark energy is uniform throughout the universe and invariable in time. This is 
the simplest explanation for dark energy, which was introduced by Einstein for building 
a static universe but has a good fit with the available data of the current cosmological 
observations. If dark energy takes this form, it is a fundamental property of the universe. 
The cosmological constant thus faces two fundamental problems in physics, namely the 
fine-tuning and coincidence puzzles. 

The late-time cosmic acceleration may alternatively be driven by a dynamic dark energy 
which could be a time evolving and spatially dependent scalar field. Lots of such dynamic 
dark energy models have been proposed, which are roughly classified into three categories: 



quintessence 



15 



-|2l|. phantom 



22- 



29- 



In quintessence models, a scalar 



and quintom 

field if with a canonical kinetic energy and a self-interaction potential energy V{ip) is sup- 
posed to be minimally coupled to Einstein gravity. In a flat Robertson- Walker background 
the quintessence scalar behaves as a perfect fluid with an evolving equation-of-state (EoS) 
parameter u = p/p lying in the range — 1 < < 1. In phantom models the quintessence 



is replaced by a ghost scalar of which the 
the no-go theorem proposed in Ref. 
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dnetic energy is negative and a; < — 1. Due to 
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40|, the model buildings of quintom dark 



energy are generally very complicated 36|]. The simplest quintom model is composed of 



two scalar fields, one is a quintessence scalar and another a phantom |29l . |30|. Quintom 
models characterize themselves by the property that the effective EoS parameter can cross 
the cosmological constant boundary u = —1, which makes them to fit the observational data 
better [36 1. 



Crossing the a; = — 1 divide in a dynamic dark energy model is bewitching. However, the 



emergence of a phantom mode with negative kinetic energy in quintom models brings about 
great embarrassment in understanding it. The consistence coming from the Null Energy 
Condition in physics requires the kinetic energy of a normal scalar field not to be negative, 
otherwise the theory might be unstable and unbounded. Therefore, it is worthwhile to study 
the mechanism of removing the phantom field from the quintom models. In fact, there has 
lots of such attempts to investigate the possibility of u = —1 crossing in quintessence like 



models 
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lN43| . It has been empirically realized that to cross the u = —1 barrier 
and remove ghost mode at the same time, the model building should be involved in either 
modifying the general theory of Einstein's relativity or introducing some higher derivative 
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43| of a 



terms for the scalar fields. For example. In the so-called Galileon cosmology 
scalar field, the higher derivatives of operators are introduced into the Lagrangian but the 
equation of motion of the scalar remains of the second order. The Galileon models can have 
oj = —1 crossing without ghost modes involved. It goes without saying, however, that these 
models are generally very complicated to deal with. 

Of the modification attempts beyond Einstein's general theory of relativity, there is an 



extended theory of gravity (ETG) based on the so-called Lyra's geometry |44j-|48l| . As is well 
known, the general theory of relativity is a theory of gravity built on (pseudo-) Riemannian 
geometry. Lyra's geometry is a modification of Riemannian geometry where a gauge function 
is introduced. Due to the presence of this gauge function on the structure-less manifold, 
an extra geometrical ingredient, i.e., the displacement vector arises in Lyra's geometry. 
It is remarkable that the connection in both Riemannian and Lyra's geometries are metric 
preserving. The extended theory of gravity in Lyra's geometry is much motivated by the 
fact that it could predict the same effects as Einstein's general relativity within observations 



limits, as far as the c 



the field equations [48 



assical Solar System, as well as tests based on the linearised form of 



59l |. The extended theory of gravity on Lyra's geometry distinguishes 



itself by the fact that it is a scalar-tensor theory of gravity, where the scalar field is not alien, 
but intrinsic to the geometry Moreover, in the so-called normal gauge |60l-l64| . a 



constant displacement vector can play the role of a positive cosmological constant (in the 



presence of other cosmic matter ingredients) 65l-l67l|. which is in contrast to general relativity 



where the cosmological constant must be added in an ad hoc manner into the gravitational 
field equations. 

In this paper, we study the a; = — 1 crossing possibility in some quintessence models in 



the framework of the ETG in Lyra's geometry. Despite the impossibihty for a canonical 
quintessence model to cross the phantom divide u = —1 in Einstein's gravity in pseudo- 



Riemannian geometry 
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371 ]. the existence of a displacement field in the ETG in 
Lyra's geometry does probably modify the effective distribution of the cosmic fluids so 
that the EoS parameter of the quintessence scalar may cross this boundary. Aimed at 
finding the exact solutions of the modified Fridemann equations in a fiat Robert son- Walker 
background, we propose several candidate interactions between the quintessence scalar and 
the displacement field which have simple mathematical expressions. For some of these 
possible interactions, crossing u = —1 barrier for quintessence models in Lyra's geometry 
is available. Unfortunately, crossing this phantom divide in these models will, without any 
exception, give rise to the instability of the relevant perturbations. The reliable quintessence 
scenarios in Lyra's geometry are still those in which —1 < u < 1. 

The paper is organized as follows. Section II begins with a brief introduction to mod- 
ified Einstein equations in the considered ETG and its application to the cosmology of 
a quintessence model in a fiat Robertson- Walker background spacetime. By analyzing the 
equation of motion of the quintessence scalar, we determine phenomenologically several can- 
didate interactions between this scalar and the displacement field. In Section III we study 
the quintessence cosmology for each of the candidate interaction terms. The self-interaction 
potential of the quintessence scalar is not given a prior, which is defined during the process 
solving the modified Friedmann equations, motivated by the requirement to have closed form 
solutions to these equations. Among the four quintessence models proposed, three of them 
naively allow = — 1 crossing. By requiring the squared sound speed of the quintessence 
scalar preserves finite and non-negative during its evolution, the possibility for u crossing 
the phantom divide is excluded. It turns out that the reliable quintessence models in Lyra's 
geometry are also characterized by inequalities — 1 < w < 1, similar to those in pseudo- 
Riemannian geometry. We conclude in Section IV with a summary of the results and some 
remarks. For simplicity we work in the Planck units c = h = = 1 throughout the paper. 

QUINTESSENCE AND ACCELERATED EXPANSION 

The quintessence models in an extended theory of gravity (ETG) in Lyra's geometry, 
in the so-called normal gauge, is described by the following modified Einstein gravitational 
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field equations |47l |: 



where, 

3 . 1 . 



= -2 - 2^?/^-/5a/3"J , (2) 

is an intrinsic geometrical stress tensor, corresponding to the existence of the displacement 

vector which emerges from the integrability condition of length of a vector under parallel 

transport [i^ ^ • T^,^ is the stress tensor of the quintessence scalar ip which is assumed to 

have a self-interaction potential V{ip) and be canonically coupled to gravity, 

1 
2 



= - -g^uVxfV\ - g^uV{<f) (3) 



We further assume that, at present epoch, the quintessence scalar cp dominates over other 
cosmic fluids such as baryonic dust and radiation. The displacement vector is allowed to 



be a time-like 4- vector field 
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62| 



= 0, 0, 0) (4) 



Its unique non- vanishing component P{t) can either be a constant or time- dependent. 

In a fiat Friedmann-Lemaitre-Robertson-Walker background ds'^ = —dt^ +a^dx^, the 
modified Einstein equations given in Eq.([T]) become: 

2,H' = \iP' + V{^)-\m (5) 

2H + ?>H' = -]^^' + V{v) + \m (6) 

where an overdot denotes derivative with respect to the cosmic time t. The equation of 
motion of the quintessence scalar, which comes from the Bianchi identities of modified 
Einstein equations, i.e., from the compatibility of Eq.Q with Eq.([6]), reads, 

(y3 + 3i/0 + V^) = ^(e + QHO) , (7) 

In Eq.([7]) V^ip = ^ and 6 = (3'^ [0 is also referred to as the displacement field). Eq.(I7j) implies 
that the displacement field and the quintessence scalar ip interact as the universe evolves. 

The mixture of the quintessence scalar ip and the displacement field 9 is conventionally 
viewed as a perfect fluid, whose energy density and pressure are defined by, 

p = \^" + y(v) - \m (8) 
p = \^"- y^v) - \m (9) 



where = i^Lp^ + V{lp) and 



V{}p) are respectively the energy density and 



pressure of the quintessence scalar 99. The real displacement vector (or positive ff) that 
is a necessary geometrical ingredient in Lyra's geometry, nevertheless, behaves as an exotic 
cosmic matter with negative energy density and negative pressure, pe = Ve = ~|^- the 
absence of quintessence scalar, the effective EoS parameter of the displacement field is equal 
to We = 1, which corresponds to the so-called stiff fluid 66|. When the quintessence scalar 



exists, the displacement field would probably play the role of a phantom 
mixture probably behaves as an effective quintom 29 
expansion of our universe. 

The effective EoS parameter of the mixed fluid is. 



TMm, so that the 



30| to cause the late time accelerated 



u :- 



P 



- AV{ifi) - W{t) 



(10) 



p 202 + AV{v) - Wit) 

If there were no displacement field 9 (as in general relativity based on pseudo-Riemannian 
geometry), the EoS parameter u would only include the contribution of quintessence scalar 



<y9, and — 1 < w < 1. This is not the case in the ETG in Lyra's geometry [66|, |67|. In 
a Lyra manifold with the positive displacement field 9, the universe evolves between u = 
1 (stiff matter), where either the kinetic term dominates or both kinetic term and the 
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24i\ . where u < —1, provided that 



displacement field dominate, and a phantom regime 
the kinetic term of the scalar field is negligible. The interesting possibility for u crossing the 
cosmological constant boundary seems plausible. To clarify such a possibility, the mechanism 
describing the possible interactions between the quintessence scalar and the displacement 
field is required. Owing to the lack of such a mechanism, in this paper, we determine the 



relevant interaction terms by a simple dimensional analysis. The interaction strengt 



inevitably encoded into a few phenomenological parameters. Following other works 68 
we interpret Eq.(I7]) as an effective energy conservation equation, and recast it as: 

(0 + 3Hip + V^) = 



h IS 



69|, 



(11) 
(12) 



where ^ stands for the required interaction. By dimensional analysis, ^ is generally of the 
form J2 = aipi + 0:2/52 + ■ ■ ■ , with pi the different energy components. The coupling param- 
eters Oj have the dimension of Hubble parameter H or 0, and should not be simultaneously 
set to zero. In the models under consideration the possible candidates of pi are ifP', 



and 9. We will simply assume that, 

36 Q~ 

^ = 3ci/02 ^ ^ ^ QfHViif) (13) 

where c, 6, 6 and / are some dimensionless coupling constants. The choice for the interaction 
terms in Eq.( IT3|) allows us, not only to explain the present accelerated expansion of our 
universe, but to have closed form solutions of the modified Friedmann equations also. 

We now proceed to study the cosmological consequences emerging from each of these 
interactions. Instead of t, we will use the e-folding number x = In a as the time variable for 
convenience, with x = representing the present time (a(0) = 1). The modified Friedmann 
equations in Eqs.([5]) and ([6]) are recast as: 



dx 



+ 6H' = 2V (14) 



o-lfHZ) ^tv-^f' (15) 



2 

3 \dx J 3 

Similarly, Eqs. ffTT]) and f|T2|) that can be viewed as the equations of motion of quintessence 
scalar and displacement field become: 



dip 
dx 



dx"^ \ 2 dx J dx 



J (16) 



4:\dx 

where. 



Relying on the fact that Eq. ([7]) is the Bianchi identity of the modified Friedmann equations 
([5]) and ([6]), any non-degenerate combination of three of equations f|T^ . f|T5|) . f|T6|) and f lT7|) 
will be mathematically equivalent. 



Case 1 c^0but6 = 6 = / = 



In this case, Eq. (ITB]l pretends to decouple from the displacement field 6'and becomes an 
effective equation of motion of quintessence scalar p, 



^2^ 

dx"^ 



+ 



3(1 - c)H' 



IdH'^'ld'^ y 
2 dx dx '"^ 







(19) 
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Equivalently, 



H 



2 1 o ^V' 







(20) 



dx"^ dx J ' \ '"^ ' dx 
when Eq. ffT^ is taken into account. In view of the mathematical structure of Eq. fl20|) . we 
assume that the self-interaction quintessence potential V{ip) is defined by, 

dip 



dx 







(21) 



in the models under consideration. The evolution of quintessence scalar in terms of efolding 
time X turns out to be independent of the Hubble parameter, 

d'^tp dip 



dx"^ dx 



(22) 



By assigning the initial conditions ip{0) = p and dip{0)/dx = q, with p, q two outstanding 
constants, we can solve Eq.( l22|) analytically. The solution reads. 



<^(x) =p + q{e^'^ - 1) 



(23) 



The constant p does not matter whether the scalar field is a quintessence or a phantom, 
because it is irrelevant to the kinetic energy of the field. So, we set p = for simplicity from 
now on. Substitution of Eq. (l23|) into Eq. fl2Tl) leads to: 



V{x) = tHq exp 



3 2 / 

-2^(1 (e 



6cx 



1) 



(24) 



In Eq. fl2^ we have used a constant r labeling the present value rH^ of the quintessence 
potential, with Hq the Hubble parameter at present epoch. The potential can be expressed 
as a closed form function of the quintessence scalar ip itself, 

3 



V{ip) = tHq exp 



-Cip{ip + 2q) 



(25) 



If c < and g = 0, such a V{ip) could be interpreted as the 



70 



;achyon potential describing the 



71| . with tHq the brane tension 



excitation of massive scalar fields on the anti-D branes 
and —3c the mass squared of the field ip in the Planckian units. However, g = is forbidden 
in our case. The physics behind the proposed potential is an open issue. 

With Eq. (!25l) . we can obtain the evolution of Hubble parameter H by solving Eq.([T] 
The result is, 

,2 



H'ix) 



Hie''''' 



r 

1 + 7^ 



3c V 30(7^ 



n 1 ^ 

c' 2 



(26) 



where r(r, 2;) is the upper incomplete Gamma function 72 1 



T(r.z) 



(27) 



Furthermore, employment of Eqs. (!25|) and f l26|) in Eq.( !T5|) yields, 



9{x) = -tHq exp 
3 



+ 2/7o (3c'g'e^^^ - 2)e- 



-6x 



+ |^.(3cVe--2)y^ 



r 



\c 2 



-r 



Vc' 2 



exp ( -cq — 6x 



(28) 



The cosmology of the mixed fluid, at the background level, is determined by its EoS 
parameter u (See Eq.l TTII]) ). For the models under consideration, 

2c r exp [6x — |cg^(e^'^^ — 11 



u{x) = 1 



3c + r 



□3cg2/2 



rQ,¥)-r(i,^e6-) 



(29) 



Although the evolution of EoS parameter u depends upon three parameters, i.e., c, q and 
r, its value uq at the present epoch is completely given by the dimensionless parameter r 
which represents the present-epoch value of the quintessence self-interaction potential, 

2r 



1 - 



(30) 



Provided r > 2, < — 1/3, the late-time accelerated expansion occurs. In particular, the 
mixed fluid in ETG of Lyra's geometry will respectively mimic the quintessence, cosmological 
constant and phantom in the Einstein's general relativity if the parameter r takes its value 
in the regions 2<r<3, r = 3 and r > 3. However, for r > 3, the EoS parameter u given 
in Eq. fl29]) increases monotonically near a; = for real parameters c and q, as seen from the 
asymptotic expansion of Eq.( 12^ at small x, 



2r 2r, , , 
1- — + — (9cV + 2r-6)a; 



(31) 



Such a u conflicts with our common sense about the universe evolution. We choose to 
abandon this possibility. The parameter r is consequently restricted within the region 2 < 
r < 3 for these models, in turn the aspired w = — 1 crossing is unavailable. 

It is interesting to study the dependence of the evolution of u upon the magnitude of 
the dimensionless coupling parameter c. To this end we plot Eq.( l29p in FIG. 1. for three 



different choices of coupling c, i.e., c = 2, c = 1.5 and c = 1, and the parameters q and 
r are fixed at g = 0.05 and r = 2.75, respectively. It is manifest that the evolution of the 
universe from matter dominant era to the present acceleration phase depends weakly upon 
what the coupling constant c is. However, the value of c will strongly influence the universe 
evolution in the future. The larger the value of c is, the earlier will the universe exit from 
its accelerated expansion. 

^ 1.0 




FIG. 1. Evolution oiuj versus x for different coupling parameter c for interaction term J2 = ScHip^. 
Here we take q = 0.05 and r = 2.75 (So ~ —0.83). The solid, dashed and dotted curves 
correspond to c = 2, c = 1.5 and c = 1 respectively. 



The above is the cosmological implications of the background dynamics of the proposed 
model in Lyra's geometry. The concordance cosmology is a science based on precise ob- 
servations of which lots are tightly connected to the growth of perturbations. Thus we 
must examine the stability issue of the perturbation modes in the model under considera- 
tion. According to the linear perturbation theory, the stability of the linear perturbation 



modes during their evolution requires j > for each component fluid [3l|, l69|, l73|], where 
c^j = dpi/dpi is its squared sound speed at the background level. In Lyra's geometry, 
the squared sound speed of the displacement field is always definitely positive. In fact, 
cj. g = ug = 1. On the other hand, 

2rexp [6x — \cq^{e^^^ 



- 1) {3 + ^^e-«V. [r(i, \cq-) - r(l, fcg^e-^)] } 
10 



The squared sound speed c^^ varies continuously with respect to e-folding time unless c = 1. 
For a positive c (c 7^ 1), cj.^^ — )■ 1 when x — j- ±00, while, 



^ 2r [1 + (6 - 2r - 9cV)x] 



at X ~ 0. Recall that 2 < r < 3, the squared sound speed c^^ of the quintessence scalar 
diverges or takes negative values for < c < 1, the corresponding perturbation modes are 
violently unstable and do not have any physical significance. The stability of perturbation 
modes is also sensitive to the initial velocity q of the quintessence scalar. To ensure a finite 
and non-negative Cg the coupling constant c should be restricted to the region c > 1, and 
at the same time the coefficient of x in the RHS of Eq. fl33p should be set to zero. Therefore, 
the model is physical acceptable only if c > 1, g = ±^a/6 — 2r and 2 < r < 3. 



Case 2 6^0butc = 6 = / = 0: 



In this case, Eq.l TTB]) is translated into. 



dx^ ^/Q \ dx J 



llJb I O 







(34) 



Similar to Case 1, we further assume that the self-interaction quintessence potential V^(v5) 
satisfies the constraint condition, 



(35) 



Consequently, the evolution of quintessence scalar in the models under consideration is also 
fictitiously independent of the evolution of Hubble parameter. 

The solution of Eq. (l36|) which satisfies the initial conditions (y9(0) = and dLp{Qi)/dx = \/6q 
reads, 

(p{x) = — — In [cosh(6x) — gsinh(6x)] (37) 
Plugging Eq.([37D into Eq.([35]) yields, 

6(1 — g^) sinh(6a;) 



V{x) 



[cosh(6a;) — gsinh(6a;)]' 



■ exp 



— 6x + 



cosh(6x) — gsinh(6a;) 



(38) 
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where r is an integration constant which is, as before, used to specify the present-epoch value 
of quintessence potential, ^(0) = rH^. Different from Case 1, for the current models, it is 
difficult to express the quintessence potential as a closed form function V{ip). Fortunately, 
this does not affect our investigation to cosmology. With (1551) . we can easily solve Eq. flT^ 
and obtain the evolution of Hubble parameter in these models. 



H\x) = 
and then, 



+ 



3(1 -g^) 3(1 -g 



exp 



6(1 -g2 



sinh(6x) 



cosh(6a;) — gsinh(6a;) 



9{x) 



4i72e-6^(3g2 + r-3) 
3 [ cosh(6x) — q sinh(6x)] ' 



(39) 



(40) 



Therefore, 



u = 1 



2(1 — g^) rexp 



6(l-g^) 



sinh{6x') 



cosh{bx)—q smh{bx) 



[cosh(6a;) — gsinh(6a;)]' 



3(1 — g^) — r + r exp 



6(l-g^) 



smh{bx) 



cosh(fex)— g sinh(fex) 



(41) 



The evolution of the effective EoS parameter depends upon three parameters, i.e., b, q and 
r. However, as in Case 1, its present value uq is only relevant to r. 



2r 

1 

3 



(42) 



Provided r > 2, uq < —1/3, the late-time accelerated expansion occurs. Notice that the 
asymptotic expansion of Eq.( HT]) at small x is. 



1 



2r 



Arx 



{l-q^)-\{r-bq) 



(43) 



So long as r < 6g + 3(1 — g^), the EoS parameter uj will decrease monotonically near 
X = for real parameters b and g, implying that the w = —1 crossing is possible. For a 
given coupling constant 6, the initial "velocity" g of the quintessence scalar have to take its 
value in the region (6 — -\/P~+T2)/6 < g < (6 + ^/h^ + 12)/6 to guarantee the inequality 
2 < r < fog + 3(1 — g^). In FIG. 2. we plot the evolution of EoS parameter u in Eq.f l4ip for 
three different choices of coupling 6, i.e., b = 2.6, 6 = 3 and b = 3.6, and the parameters g 
and r are fixed at g = 1 and r = 2.05, respectively. The remarkable difference between the 
present case and Case 1 is that the EoS in the present models can cross the phantom divide 
a; = — 1. Is this quintom scenario reliable? 
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1,5 r 




-1.0 - \ 



-1.5 L 

FIG. 2. Evolution of co versus x for different coupling parameter b for interaction term ^ = 
SbOif /2V6. Here we take q = 1 and r = 2.05. The solid, dashed and dotted curves correspond to 
b = 2.6, 6 = 3 and b = 3.6 respectively. 



It has been pointed out 



31 



36| that a viable quintom scenario can not be realized only by 



virtue of the parameterization of EoS. The stability of the relevant perturbation modes must 
be ensured also. In other words, we have to guarantee c^^ > at background level. For 
simplicity we only consider a special case q = 1- In this case, the solution to the background 
dynamics reduces to: 



V{x) 



X 



1 + -(e 
b 



2bx 



-V. 



[X 



-rij2e2{''-3)- 



2bre 



2bx 



3[6 + r(e2''^ - 1) 
The energy density and pressure of the quintessence scalar read. 



(44) 
(45) 
(46) 

(47) 
(48) 



Hence, 



tt2 
b 



[3(6 - r) + (3 + b)re^'"'] 



HI 



e-^"[3(6-r) + (3-6)re2''"] (49) 



dp^ _ (3 - 6)Ve2''^ + 9(6 - r) 
dp^ ~ (9 - 62)re2^^ + 9(6 - r) 
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(50) 



The physical acceptance requires w < — 1/3 for x > and c^^ > 0. Obviously, both 
inequality can be satisfied if2<r<6<3. With respect to such a parameter constraint, 
however, 1 > > 1 — 26/3, the quintom scenario where the EoS of mixed fiuid can cross the 
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cosmological constant boundary is forbidden. It seems that the no-go theorem 
40| is valid also for a generic q so that a reasonable quintom scenario is unavailable in the 
present model. 




FIG. 3. Evolution of Cg versus x for different coupling parameter h for interaction term £2 



36610/2^/6. Here we take g = 1 and r = 2.05. The solid, dashed and dotted curves correspond to 
6 = 2.6, 6 = 3 and 6 = 3.6 respectively. In the first two cases, is finite and positive. In the 



last case, ,^ diverges during its evolution, implying the instability of the relevant perturbation 
modes. 



Case 3 6^0butc = 6 = / = 



In this case the interaction between the quintessence scalar and displacement field is ^ = 
^bH9. We choose Eqs. ffT^ . (|T6|l and f|T71) to form the set of independent equations. The 
latter two can be recast as, 



dx 



+ 6(1 



+ 3 



1 dH^\d^ 



b)9 



2H'^ dx J dx 
- 



+ 



dV 

dx 



9~ 

-be 

2 







(51) 
(52) 
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among which, Eq.( l52l) is easily to solve. By assigning the initial condition 6{0) = 4:[b—l)sH^, 
with s a dimensionless constant, we have the solution of Eq.( l52|) as follow. 



9{x) = 4(6 - l)s i/oV^^-')" 



(53) 



We further assume that in the models under consideration the self-interaction potential of 
the quintessence scalar possesses property. 



dx 2 



Under this assumption, Eq. (l^ is reduced to. 



+ 3 



1 dH^\ _^ 



dx"^ \ 2iJ2 (]^x J dx 
Substitution of Eq.([53l) into (IMI) yields, 

V{x) = Hi \r + 3fe(e^(^-i)" - 1) 



(54) 



(55) 



(56) 



where r is an integration constant. As in the previous cases, this parameter characterizes 
the present-epoch value of the quintessence potential. With Eq. fl56l) . we can obtain the 
evolution of the Hubble parameter by solving Eq. (fT^ . The result is. 



Hi 



hs 



s e 



-6x 



+ se^(^-i)^ - bs 



-6x 



(57) 



Plugging Eq.f lFTI) into gives, 

d(p 
dx 



' (58) 

'3+{r- 3fe)(e6^' - 1) + 3s(e6^^- - 1) 

where the integration constant q is related to the initial "velocity" of the quintessence scalar 
by d(^{0)/dx = q/y/S. The cosmology of the considered models at the background level is 
described by the following effective EoS parameter: 

r + 36s(l - e6(^-i)^) 



CO 



1 - 2e 



(59) 



3 + (36s - r)(l - e6^0 - 3s(l - e^^^) 

The EoS parameter seems not to depend upon the choice of the initial "velocity" q of the 
quintessence scalar, but upon the initial value s of the displacement field instead. This 
is, however, merely an optical illusion. The consistence of the above results with Eq. ffTSj) 
requires. 



{b-l)s 



g2 + 6(r - 3) 
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(60) 
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The present-epoch value of EoS parameter is still given by the same formula as either Eq. (13T]) 
or fH2|) . which depends only upon the parameter r, 

c^o = 1 - f (61) 

Provided r > 2, the universe is destined to enter a late-time acceleration phase. 

To examine the stability of the relevant perturbation modes, we calculate the squared 
sound speed of the quintessence scalar. The result is. 



., 3 - r + 3 6 - I s + 3b(b - 1 se*^*^ , , 

Ci^= ^ (62) 

3-r + 3(6- l)s-36(6- l)se6*^ 



As long as (6 — l)s ^ 0, cj.^^ either diverges or becomes negative during its evolution, 
which will give rise to the unstable and then the unacceptable perturbation modes. On 
the other hand, when (6 — l)s = 0, which occurs for either 6 = 1 or s = 0, c^^ = 1, the 
relevant perturbation modes might evolve stably. For 6 = 1, the solution to the background 
dynamics reduces to: 



<p{x) = -f-t.nh^' f ;_+;^" (63) 

V{x) = rHl (64) 

= -Hi [r + (3 - r)e-H (65) 
3 - 

e{x) = (66) 

and in particular, 

2(3 -r 

3 — r + re' 

In this case, the displacement field is effectively absent, and the self-interaction potential of 
the quintessence scalar plays the role of the cosmo logical constant. The potential parameter 
r must be restricted to the region 2 < r < 3, otherwise u will diverge during its evolution. 
Consequently, crossing the phantom divide u = —1 in the present model is practically 
impossible. Fig. 4. shows the evolution of EoS parameter in Eq.( l671) for three choices of 
the potential parameter r, i.e., r = 2.6, r = 2.9 and r = 3, with the coupling constant 
fixed at 6 = 1. It appears that the accelerated expansion in these models will last for a very 
long time. Besides, the larger the parameter r is, the more closely the quintessence scalar 
resembles the cosmological constant. 
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X 



FIG. 4. Evolution of oj versus x for different potential parameter r if the interaction term is 
J2 = ^bH6 with 6=1. The solid, dashed and dotted curves correspond to r = 2.6, r = 2.9 and 
r = 3, respectively. 



Case 4 / ^Obut6 = 6 = c = 0: 



In this case, the interaction between the quintessence scalar and the displacement field is 
assumed to be proportional to the quintessence self- interaction potential, ^ = 6fHV{ip). 
Under such an assumption, Eqs.(lT6l) and ( |T7I) become. 



H 



.(Pip 



dip 
dx 

^ + 69-8fV = 
dx 



, . ldH'^\dip 



2 dx J dx 



+ ^-6/^^ = 

dx 



As before, we further assume that the quintessence potential is defined by condition, 

dx 



(68) 
(69) 

(70) 



This implies that the quintessence self-interaction potential in the considered models is of 
the form, 

V{x) = rHy^'' (71) 

where a real parameter r is used to characterize the present-epoch value of the potential 
and Hq stands for the present value of the Hubble parameter. Substitution of Eq. (!7T!) into 
Eq.(fTiD yields. 



2-6x 



^3(/ + l)^'' 



(72) 
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Furthermore, we can obtain the evolution of the quintessence scalar by plugging Eqs.(17T 
and ([72]) into Eq.(|68l), 



V3(/ + l)(3 + 3/- 



tanh-S/ /^^ + ^^ -tanh-^ 
3 + 3/ — r 



'3(/ + l)+r(e6(/+i)--l) 



3 + 3/ - r 



(73)" 



where the parameter q is the integration constant which can be interpreted as the initial 
velocity of the quintessence scalar, dip{0)/dx = q. Finally, the evolution of the displacement 
field is obtained from Eqs.([I5]), (I7T]). (172]) and (175]). 



0{x) = ^Hle-'^ 



q' -6 + 2r 



Jg6(/+l)x' _|_ I 

7+^ , 



(74) 



The effective EoS parameter of the quintessence scalar and the displacement field in the 
models under consideration reads, 

2(/ + l)re6(/+i)- 



1 



(75) 



3(/ + l) + r(e6(/+i)--l) 

The present-epoch EoS parameter takes the same formula as those in the previous three 
cases, 

2r 

coo = l-- (76) 



So the late-time cosmological acceleration is available in these models if r > 2. Different 
from the EoS parameter in the present case which depends only upon the coupling constant 
/ and the potential parameter r, the squared sound speed of the quintessence scalar depends 
also upon the initial velocity q of the quintessence scalar, 

.2 , 2/re'^(^+^)'' 



(77) 



q2 _ 2/re6(/+i)^ 

Stability condition ^ > requires / = and q ^ 0. Therefore, the model is physically ac- 
ceptable only if there is no interaction between the quintessence scalar and the displacement 
field. When / = 0, cj.^^ = 1, the solution to the background dynamics reduces to: 



dx ~ y/3 

Vix] 



VSq 



r + re 



6x 



rHl 



e{x) 



1 



r + (3 — r)e 



—6x 



^oe"'"(g'-2r-6) 



(78) 
(79) 
(80) 

(81) 



18 



Because of 6{x) > 0, the value of parameter q is restricted to either q > \/2{3 — r) or 
q < -^2{3-r). For / = 0, Eq.([75D reduces to: 



00 = 1 



3 — r + re' 



6x 



(82) 



To have w < — 1/3 at x = 0, r > 2. To avoid the possible divergence of u during its 
evolution, r < 3. FIG. 5. illustrates the dependence of the EoS parameter in Eq.fl82l) upon 
the potential parameter r. Obviously, — 1 < w < 1, crossing the phantom divide a; = — 1 is 
prohibited once more in the quintessence model in Lyra's geometry. 




FIG. 5. Evolution of a; versus x for different potential parameter r if there is no interaction between 
the quintessence scalar and the displacement field. The solid, dashed and dotted curves correspond 
to r = 2.1, r = 2.9 and r = 3, respectively. 



CONCLUSION 

In this work, we have established four classes of the quintessence models in ETG in 
Lyra's geometry. The classification of these models depends upon how the quintessence 
scalar interacts with the geometrical displacement field 6 (or The mixture of in- 

teracting quintessence scalar and the displacement field supplies as a cosmological perfect 
fluid which can cause late time accelerated expansion of our universe. Owing to the subtle 
choices of the quintessence self-interaction potential, all quantities relevant to the study of 
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cosmology, including the Hubble parameter, the displacement field, the time derivatives of 
the quintessence scalar and the potential itself, are expressed as closed form functions of 
efolding time x = ln(a), and so is the effective EoS parameter ui of the mixed fiuid. The 
evolution of oo is different for the quintessence models of different classes, which depends also 
upon the present values of the time derivative of quintessence scalar, its potential and what 
the coupling constant is. However, today's u) does only depend upon the present-epoch value 
of the quintessence potential. The appearance of the displacement field in Lyra's geometry 
improves greatly the late-time evolution of the quintessence-displacement mixture, however, 
crossing the phantom divide a; = — 1 is still forbidden in these models by the necessary 
condition ^ > to ensure the stability of the cosmological perturbations. Establishing a 
reliable quintom scenario remains a challenge in ETG in Lrya's geometry . 

In the proposed quintessence models in ETG in Lyra's geometry, we have defined the 
quintessence scalar by making some careful choices for its self-interaction potential and 
interaction terms with the geometrical displacement. These choices, in this paper, are mainly 
motivated by the requirement to obtain the analytical solutions of the modified Priedmann 
equations. More important issue that remains unsolved is to investigate the physics behind 
these choices. Why the EoS parameter of the quintessence scalar at present epoch depends 
only upon the parameter r for all four kinds of models is also a mystery. In addition, the 
characteristic behaviour that the effective EoS parameter u decrease from a; 1 in the past 
confiicts with the well-known Big-Bang diagram. In our a; ~ x figures the dust-dominant 
phase does not form an expected plateau. To be more realistic, the pressureless cold dark 
matter components have to be introduced into the model buildings also. 
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